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We construct a universal set of high fidelity quantum gates to be used on a sparse bipartite 
lattice with always-on Ising couplings. The gates are based on dynamical decoupling sequences 
using shaped pulses, they protect against low-frequency phase noise, and can be run in parallel 
on non-neighboring qubits. This makes them suitable for implementing quantum error correction 
with low-density parity check codes like the surface codes and their finite-rate generalizations. We 
illustrate the construction by simulating quantum Zeno effect with the [[4, 2, 2]] toric code on a spin 
chain. 



Quantum error correction (QEC) makes it theoreti- 
cally possible to perform large quantum computations 
with a finite per-qubit error rate[TH3]. In practice QEC is 
extremely difficult since the corresponding error probabil- 
ity threshold is small[lH5]- When only local interactions 
between the qubits are allowed, the estimated threshold 
value is the highest, around 1%, for toric and related sur- 
face codes[Slini[in]- However, as to how to implement the 
operations efficiently and with the required accuracy, is 
an open question. 

Qubits with always-on couplings are a natural model 
for several potential quantum computer (QC) architec- 
tures currently under investigation. On general grounds, 
compared to their counterparts with tunable couplings, 
such qubits can be expected to have better parameter 
stability and longer coherence times. In addition, over 
sixty years of development in nuclear magnetic resonance 
(NMR) yielded an amazing degree of control available 
to such systems [m [T^. Related techniques based on 
selective dynamical decoupling (DD) of parts of the sys- 
tem Hamiltonian with carefully designed pulse sequences 
have been further developed in application to quantum 
computation fTMTO] . 

While NMR quantum computation is not scalable [20]. 
it still holds several records for a number of coherently 
controlled qubits[12]. However, some of these have been 
achieved with the help of the strongly-modulated pulses, 
computer-generated single- and multi-qubit gates tai- 
lored for a particular system Hamiltonian piH^ . While 
such gates can be used in other QC architectures [25] . 
they may violate the scalability. 

On the other hand, NMR-inspired techniques like DD 
can also be used to control large systems with local inter- 
actions, with the scalability achieved by designing pulses 
and sequences to a given order in the Magnus series [26] 
on small qubit clusters [?fl 128]. DD is also excellent in 
producing accurate control for systems where not all in- 
teractions are known as one can decouple interactions 
with the given symmetry f 291 150]. Moreover, DD works 
best against errors coming from low-frequency bath de- 
grees of freedom which tend to dominate the decoherence 
rates, and it does not require additional qubits. In short. 



DD is an excellent choice for the first level of coherence 
protection; it's use could greatly reduce the required rep- 
etition rate of the QEC cycle. 

This is well recognized in the community and applica- 
tions of DD for QC are actively investigated by a number 
of groups. However, most publications on the subject il- 
lustrate general principles using a single qubit as an ex- 
ample, leaving out the issues of design and simulation 
of scalable approaches to multi-qubit dynamical decou- 
pling. While the techniques for larger systems exist, they 
typically yield longer decoupling sequences [TT ] [29 ] l3T] . 

The goal of this work is to provide a scalable 
benchmark implementation of a universal set of accurate 
gates using soft pulses for a system with always-on qubit 
couplings. Specifically, we construct one- and two-qubit 
gates with built-in DD-protection against low-frequency 
phase noise for a sparse bipartite lattice of qubits with the 
nearest-neighbor (n.n.) Ising couplings. The constructed 
gates use finite-amplitude shaped pulses which can be 
implemented experimentally. They are scalable, in the 
sense that the same construction works for an arbitrary 
lattice, and they can also be executed in parallel for dif- 
ferent qubits and/or qubit pairs. This makes them ideal 
for implementing QEC with quantum low-density parity 
check (LDPC) codes[32l [33], in particular, the surface 
codes and their finite-rate generalizations [5] 134] l35]. In 
the limit of very slow (classical) bath the gates are ac- 
curate to second order in the Magnus expansion, mean- 
ing that their infidelity scales as sixth or higher powers 
of the coupling, in units of inverse pulse duration. We 
demonstrate the accuracy of the constructed gates by 
simulating the quantum Zeno effect [5S1 157] repeatedly for 
the [[4, 2, 2]] error-detecting toric code on an Ising chain. 
The simulations are done with five qubits, using classical 
correlated noise as a source of dephasing. 

Two techniques are essential to our work. First, the 
use of NMR-style self-refocusing pulses [?ni551fln] . which 
(to a given order) work as drop-in replacement for hard, 
(5-function-like pulses. In simulations, we use the second- 
order pulses designed and characterized in Refs. [27l [28l 
139] andim The second technique is the Eulerian path 
construction for generating accurate DD sequences [TB], 
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and its extension, the dynamically corrected gates [5T1H^ 
|43] which allow for the construction of composite pulses 
accurate to a given order of the Magnus expansion. 

We construct our gates for a collection of qubits 
arranged on an arbitrary sparse bipartite graph Q with 
edge set f , with an Ising coupling for every edge, 

arbitrary (within the bandwidth) single-qubit control, 

^c^^E E ^'v^'mW, (2) 

in the presence of low-frequency phase noise 

^N^^E<^» + ^B, (3) 
i 

where the bath coupling operators Bi can result, e.g., 
from low-frequency phonons, or nuclear spins, and their 
dynamics is governed by the bath Hamiltonian inde- 
pendent from a^. 

Decoherence resulting from higher-frequency bath 
modes, e.g., described by the Lindblad equation [44]. can 
be also introduced, but at later design stages, since DD 
is not effective against such decoherence. While we do 
not consider Markovian decoherence here, we mention in 
passing that the main effects of DD are the suppression 
of equilibrium population asymmetries (qubits are con- 
stantly flipped), and, with soft-pulse DD, the redistribu- 
tion of decoherence rates between the channels 01] ■ For 
example, even if dephasing is dominant for non-driven 
qubits, any sequence of finite-width pulses creates some 
longitudinal relaxation (compensated by a reduction of 
the dephasing rate). 

To construct the CNOT gate, we use the identity [JS] 

CNOTed - e'^^/^Y^X^Y^XdYde-'^/^^Yd, (4) 

where Xi = exp(— i^crf), Yi = cxp (—i|tTf), are it/2 
unitaries, and Xi, Yi denote the corresponding conju- 
gate gates [in simulations we use the equivalent form with 
YcXcYc = Zc]- 

To IMPLEMENT THE TWO-QUBIT zz-rotation gate, 
g-i7r/4cr^(Tj^ we run two period-16rp decoupling sequences 
on the sublattices A and B, VA{t) and Vsit) in Fig. [l] 
where each pulse of duration Tp is a symmetric tt pulse 
applied in the x direction. When the pulses are second- 
order self-refocusing pulses [e.g., Qi{tt) from Refs.BTland 
[Mlshown], these sequences suppress the effect of the Ising 
couplings ffs and the noise Hn to second order in the 
Magnus expansion, meaning that the effective Hamilto- 
nian is just TJbi with the error scaling as cx . This gives 
the error in the unitary matrix scaling as cx r^, and the 
corresponding infidelity scaling as cx: t^. 




time ( in units of r ) 

FIG. 1. (color online) Pulse sequences used to implement two- 
qubit z^-rotations on a bipartite lattice with Ising couplings. 
Global sequences of vr-pulses in x direction VA(t) and Vsit) 
are executed on the idle qubits of the two sublattices. These 
decouple the inter-qubit Ising couplings and also the single- 
qubit low- frequency phase noise terms. In order to couple 
two neighboring qubits, on these qubits the sequences Va, 
Vb are replaced by Vc, Vd, respectively. This produces an 
effective Ising Hamiltonian with half of the original coupling. 
Shown are Q-^iji) second-order self-refocusing pulse shapes |27l 
I39| used in the simulations. 

To turn on the coupling between two neighboring 
qubits c, d, the two sublattice sequences on these qubits 
are replaced with Vc{t) and Vd(<), respectively, see 
Fig. [T] These sequences are chosen so that the Ising cou- 
pling between these qubits is only removed half of the 
time, while the coupling to other qubits continues to be 
removed. More precisely, with second-order pulses, the 
effective Hamiltonian is ifB + ( J/4)cr^tT^-|-C'(rp ). Repeat- 
ing this sequence m times gives the system evolution 

— exp(— iacr^CT^), a = 4mjTp, (5) 

with the error scaling as oc rriTp, where the term as- 
sociated with the bath evolution is suppressed. Such 
sequences can be run simultaneously on many pairs of 
qubits as long as qubits from different pairs are not mu- 
tually coupled. 

We IMPLEMENT SINGLE-QUBIT ROTATIONS with the 

leading-order DCGs[3TJ|35], using the pulse sequences in 
Fig. [2] Again, two decoupling sequences, Vi{t) and V2{t), 
are run globally on the two sublattices; additional pulses 
are inserted for the qubits to be rotated [V3(i) in Fig. [2] 
shows an implementation of the tt/2 rotation with respect 
to Y axis on a sublattice- A qubit]. 

Nominally, DCG guarantees the first-order decoupling 
with any pulses. However, in our case, the decoupling 
sequences Vi{t), i — 1,2, do not go over the complete 
single-qubit groups. Thus, unoptimized (e.g., Gaussian) 
pulses can produce unitary errors scaling linearly with 
the pulse duration Tp; one needs first-order self-refocusing 
pulses 15^ 155] to get first-order decouphng. In the case 
of the second-order pulses (e.g., Qi{(j>) US]), the remain- 
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FIG. 2. (color online) Pulse sequence used to implement 
single-qubit rotations. The sequences of n pulses in x direc- 
tion, Vi{t) and V2{t), are executed globally on the idle qubits 
of the two sublattices. A single-qubit rotation is implemented 
as a DCG[3n I42| by adding three pulse-antipulse combina- 
tions and the stretched pulse in the shaded regions; the se- 
quence V3{t) where Q\['k/2) pulses in y direction are added 
corresponds to a (7r/2)y operation on a sublattice-A qubit. 
Such gates can be executed on any set of non-neighboring 
qubits. 

ing order- Tp errors are all proportional to different com- 
mutators [Bi^Bj] and [H^^Bi], which gives second-order 
decoupling (infidelity oc t^) when the operators Bi are 
replaced by c- numbers (cf. chemical shifts in NMR). 

These accuracy predictions are confirmed in Fig. [3] 
which shows the average infidelities for a single {■k/2)y 
rotation of qubit 3 [Fig. [sj^a)] and a complete CNOT23 
gate [Fig.|3jb)] as a function of the r.m.s. chemical shift A 
(in units of t^^), obtained numerically for a four-qubit 
Ising chain. The simulations are done with a custom 
C++ program using fourth-order Runge-Kutta algorithm 
for integrating the unitary dynamics and the EigenS 
library[?5] for matrix arithmetics. We fix the value of 
Jij = J = 7r/(1677iTp) with m = 5 repetitions of the basic 
sequence [see Eqs. ^ and Fig.jl] in the CNOT gate; 
with the addition of four single-qubit DCGs [see Fig. [2] 
the CNOT duration is icNOT = 9 x 16Tp = 144Tp. For 
small A, the infidelities are dominated by the decoupling 
accuracy of the inter-qubit interactions, while they scale 
as oc A^Tp for large A, see the graphs of the correspond- 
ing slopes in the insets. 

We illustrate the performance of the designed 
gates by simulating the quantum Zeno effect 137) us- 
ing the four-qubit toric error-detecting code [TOl EZ] . We 
used zero-mean classical stationary Gaussian stochas- 
tic processes with Gaussian correlations, {Bi{t)Bj{t')) — 
the decoherence source. These are 
obtained by applying the Gaussian filter to discrete sets 
of uncorrelated random numbers drawn from the Gaus- 
sian distribution, and using the standard cubic spline in- 
terpolation with the result. 

The [[4, 2, 2]] toric code is a stabihzer code [Ml |49] 
encoding an arbitrary state of fc = 2 qubits into a 2^- 
dimensional subspace Q of the 4-qubit Hilbert space. The 
subspace Q is the common +1 eigenspace of the two sta- 
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FIG. 3. (color online) (a) The infidelity of the gate in Fig. [5] 
and (b) the infidelity of the complete C-not gate between 
qubits 2 and 3 of a four-qubit n.n. Ising chain as a function of 
r.m.s. chemical shift A (pulse sequences not shown). For small 
A the errors saturate because of the fixed interqubit coupling 
J; for large A the errors are dominated by the chemical shifts. 
The insets show the corresponding slopes; the slope of ap- 
proximately 6 [infidelity oc (rpA)''] indicate that the on-site 
chemical shifts are decoupled to quadratic order. Analogous 
calculations with first-order pulses 27, 39 Si{tt) and Si{tv/2) 
give Ist-order [infidelity oc (rpA)*] decoupling and two or- 
ders of magnitude higher infidelities (not shown); Gaussian 
pulses increase infidelities by up to five orders of magnitude 
depending on the width. 

bilizer generators, Gx = <j\a^a^a\ and — afa2cr^al. 
We use the following explicit map (up to normalization) 

|oo) = [0000) + [nil), |oi) = |ooii) + |iioo), 

|io) = [0101) + [1010), = [0110) + |iooi), (6) 

where digits with tildes indicate the states of the logical 
qubits. An application of any single-qubit error, i.e., a 
Pauli operator af, /j, = x,y,z, takes the encoded wave- 
function to one of the three orthogonal subspaces, where 
one or both of the eigenvalues of Gx , Gz (these eigenval- 
ues form the error syndrome) equal —1. The code has 
distance d = 2 since some two-qubit errors, e.g., cticI, 
act within the code space and cannot be detected. 

In the presence of the error Hamiltonian ([3]), to lead- 
ing order in the perturbation, the original wavefunc- 
tion ■00 evolves into a superposition of orthogonal terms 
^0 li^o) + ^i'^i IV'o): where, in general, the coefficients 
are operators acting on the bath degrees of freedom, and 
the state fidelity is = Tre ( AJAqPb), where trace over 
bath degrees of freedom with the density matrix pB is 
taken. With the same accuracy, F is also the probability 
that the measurement returns Gx = Gz = 1- For weak 
perturbation, and for times t small compared to the bath 
correlation time Tc, the infidelity I — F scales quadrat- 
ically with t; thus frequent projective measurements of 
the generators ensure preservation of the wavefunction 
with high probability (quantum Zeno effect [551 157] ). 

Using the constructed DD-based gates, we simu- 
lated the encoding/decoding and ancilla-based stabi- 
lizer measurement circuits (Figs. |4j [sL where we 
use standard quantum circuit notations ]49l 150) . The 
Hadamard gate was implemented using the identity 
Uh = e"*'^/^ exp(i|(Tj^) exp(i|(Tx), and a projective 
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single-qubit measurement was implemented simply as an 
instantaneous projection to the state |0) of the ancilla. 
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FIG. 4. Encoding circuit implemented on a spin chain using 
four exchange gates (implemented with three CNOTs each), 
five CNOT gates, and a Hadamard gate H. Input qubits ii, 
12 can be in an arbitrary two-qubit state, on the output the 
circuit returns an equivalent linear combination of the states 
in the code, see Eq. Q, using qubits Sj, j — 1, ... ,4, and 
an ancilla initialized for the stabilizer measurement circuit in 
Fig. [5] The decoding is done by reversing this circuit. 
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FIG. 5. Measurement circuit implemented on a spin chain 
with a shuttling ancilla. The ancilla is first shuttled up for 
the Gz measurement and then shuttled down for the Gx mea- 
surement (see text). The entire circuit is repeated for every 
Zeno cycle. 

Samples of the simulation results are shown in Figs. [6j 
[7] which show the time-dependence of average infideli- 
ties and accumulated circuit failure probability, averaged 
over 20 instances of the classical correlated noise. The 
time axis starts at the end of the encoding circuit, see 
Fig. |4] Closed symbols represent the computed infidelity 
after the syndrome measurement (two points per circuit 
in Fig. [5]) ; open symbols at the end of the final decoding 
(reverse of the circuit in Fig. |4]) . Blue squares stand for 
complete simulation, red circles have the pulses but no 
projection (dynamical decoupling is done but no Zeno ef- 
fect happens, increasing the infidelity in our simulations 
by up to an order of magnitude), and black triangles 
where neither the decoupling pulses nor ancilla measure- 
ments are applied. 

It is clear that the dynamical decoupling and the Zeno 
effect are both contributing to improving the fidelity, and 
both get better with decreasing noise. This indicates 
that the errors contributing to the infidelity of the con- 
structed gates are not dominated by high-weight errors 
which would be undetectable by the code|51). 



FIG. 6. (color online) Infidelities during the Zeno cycle for 
difi'erent noise correlation times as shown, with the same noise 
amplitude a = 10~^/rp. The different curves correspond to 
cases where no pulses are applied (NP), DD pulses are applied 
but no measurements are done (NM) , and with the syndrome 
measurements (WM). Note that the axis for the cumulative 
success probability (SP) of no error is on the right. 
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FIG. 7. (color online) Same as Fig. [g] but with fixed noise 
correlation time Tc = 128rp with the r.m.s. noise amplitudes 
a as indicated. 

The infidelity sharply increases with shorter noise cor- 
relation times; this results from the asymmetry of the 
DCGs, see Fig. [2] We have also constructed [51 sym- 
metrized DCGs which give second-order decoupling for 
arbitrary bath operators Bi when used with the pulses 
constructed in Ref. SOj we expect such gates to have a 
much better accuracy for smaller noise correlation time, 
down to the gate duration, Tc > 32rp. 

In conclusion, we implemented a universal set of 
one- and two-qubit gates for a system with always-on 
qubit coupling. The gates are based on DD techniques 
and have an added benefit of protection against low- 
frequency phase noise. 

One application of thus constructed gates is for imple- 
menting toric codes on square lattice, where one sublat- 
tice would be used for actual qubits, and the other sublat- 
tice for ancillas. This way, measurement of the entire syn- 
drome can be done in just two cycles, each of four CNOTs 
in duration, plus some single-qubit gates. The same se- 
quences would also work for an arbitrary quantum LDPC 
code, if the couplings between the qubits and the ancil- 
las form the corresponding Tanner graph |52]. In partic- 
ular, for hypergraph-product and related codes|311 |3S] 
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one can use the square lattice layout with additional 
connections [53], 
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